Effective forces induced by fluctuating interface: exact results. 
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We present exact derivations of the effective capillary wave fluctuation induced forces resulting 
from pinning of an interface between two coexisting phases at two points separated by a distance 
r. In two dimensions the fsing ferromagnet calculations based on the transfer matrix approach give 
an attractive force decaying as 1/r for large distances. In three dimensions mapping of the body- 
centered solid-on-solid model onto the 6-vertex model allows for exact solution using the bosonization 
analysis of the equivalent XXZ Heisenberg quantum chain. The exact result gives the attractive 
force which decays asymptotically as l/(rlogr). 

PACS numbers: 05.50.+q, 



Interfaces between coexisting thermodynamic phases 
exhibit large spatial fluctuations at temperatures above 
the roughening transition. The length scale of such fluc- 
tuations diverges with system size, so that the inter- 
face is "washed out" in the thermodynamic limit (un- 
less suitable external fields are applied), a result com- 
pletely at variance with classical thermodynamical ex- 
pectations. Nevertheless, such a striking phenomenon 
has recently been observed by a direct visual method in 
an ingeniously constructed experiment with a colloidal 
system [jj . Further, there are many implicit manifesta- 
tions in such diverse areas as biological systems the 
asymptotic behaviour of correlation functions 0] and vic- 
inal sections of crystals in the terrace-ledge-kink (TLK) 
model [1,01 ■ Here ledge fluctuations are examined, typ- 
ically when there is a short range interaction between 
such ledges, either attractive or repulsive, together with a 
long-ranged repulsion which comes about because neigh- 
bouring ledges cannot cross and so restrict each others' 
configuration space. Such a force is said to be of "en- 
tropic" or "Casimir" type. A crucial fact which will be 
used later in another context, is that the TLK model of 
a crystal surface can be mapped onto the 6-vertex model 
of Lieb [Io| and solved exactly, with various phase tran- 
sitions resulting, depending on the nature of the short- 
ranged force: attraction gives surface reconstruction, re- 
pulsion facet selection at a "magic" angle, with the facet 
disappearing by a roughening mechanism. In this Let- 
ter we consider a related interfacial problem, but one in 
which geometric limitation of the configuration of the in- 
terface manifests itself in a long range attraction. This 
we analyse in d = 2 and d = 3 by using exact d = 2 
Ising techniques in the former case and another 6-vertex 
mapping, together with some new results, in the latter 
case. This 6-vertex mapping, due to van Beijeren [f|, is 
different from the TLK one, but has the crucial common 



feature that it gives height conservation around closed 
loops, a necessary condition for the integrity of the sur- 
face. Further, it manifests a phase transition of roughen- 
ing type 0. 

It is well-known that chemico-physical forces can trap 
colloidal particles at the liquid-liquid or liquid- vapour in- 
terface at coexistence Q. This phenomenon is widely 
used, for example, to create the stable emulsions in sys- 
tems called Pickering emulsions 9]. Another scenario 
here is that of a bifunctional particle, say a sphere, the 
surface of which is divided at the equator into two hemi- 
spheres with differential wetting properties. Since the 
colloidal particles are much larger and more massive than 
the ones in the phase separating systems, they restrict the 
motion of the interface to which they are rigidly attracted 
strongly. As a limiting case, they can be treated as spa- 
tially fixed. The phase space restriction then implies, as 
we shall show later, that, were the massive particles to be 
released, then they would accelerate towards each other. 
The interest of such questions stems not only from the 
various practical applications [1, 11, [13] but also from 



the basic requirement to understand the nature of the 
effective forces between the colloidal particles [H, [I3 |. 
In particular, the important role of the aforementioned 
capillary fluctuations is clear; substantial progress has re- 
cently been made TH, 3] , but there remain a number of 

we calculate the incremental free energy of an 
interface in the d = 2 Ising ferromagnet resulting from 
pinning, which in this case can be achieved in various 
ways. The method, which carries over to our d=3 cal- 
culations, is to fix the interface at its extremities and 
to use this for the reference state. The pinned state is 
then a further restriction of the reference state in which 
two interior points are required to have the same inter- 
face displacement, but all such displacements are allowed 
with equal weight. The associated partition functions are 
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FIG. 1: Schematic plot of Ising lattice with interface pinned 
by its endpoints and in two interior point at the same level. 
The nearest-neighbor ferromagnetic coupling is Kin unit of 
kT. The transfer matrix acts in x direction. 



calculated using the standard transfer matrix technique, 



supplemented with the "domain wall state method" [ljj : 
the state \m) achieves localization at the position m on 
any line normal to the mean interface direction and can 
be written as 

\m) = M- 1/2 J2cnp-ik(m+ 1/2) ti(k) G f (fc)|$). (1) 

Here |$) is both the G(k) vaccum and the maximal eigen- 
vector of the transfer matrix V . The G^(k) are Fermi 
operators which "diagonalize" the transfer matrix for a 
lattice of width M normal to the transfer direction . 
Thinking of the transfer matrix as generating a discrete 
Euclidean time evolution, the G'(k) create "particles" 
with energy "f(k) given by Onsager's formula cosh7(fc) = 
cosh 2K cosh 2K* — cos/c, where K is the nearest neigh- 
bour coupling in units of fcgT, K* is the dual coupling 
given by the involution smh(2K) sinh(2X*) = 1. It is im- 
portant to note that |m) only localises the interface up to 
the bulk correlation length as is revealed by calculating 
the local magnetization [I?]]. The ratio of the partition 
function with two pinned states |m) separated by a hor- 
izontal distance 2n (see Fig.[IJ, denoted by Z x (n\N, M) 
to that for the unpinned system, denoted by Z(N, M), is 



Z x (n\N,M) 
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Z{N,M) (0\V 2N \0) 

x (m\V 2n \m){m\V N - m \0) (2) 
As M — > oo, bringing in expression for (m| gives 

Z x (n\N, oo) ^ [U(N — n, m)l U(2n, m) 
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where 



U(N,m) = — e- N ^ +imu . (4) 
Jo 27r 

Evaluating the sum over m in ^ gives 

Z x (n\N, oo) U{2N - 2n, 0) U(2n, 0) 
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(5) 



Finally, using the fact that 
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U(2N-2n,0) _ 2n7(0) 
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(6) 



we conclude that in the limit M, N — > oo the incremental 
free energy is given by 
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Here 7(0) is proportional to the inverse correlation 
length. The asymptotics for large n of this are deter- 
mined by the Laplace method 



p f x ^log v /2 7 r 7 (2)(0)r/ ao + 0(l), 



(8) 



where r = 2nao and ao is the lattice spacing, valid for 
n7 (2) (0) > 1, i.e., r > a /7 (2) (0), where l/7 (2) (0) is 
the interface stiffness. The stiffness comes in because the 
surface tension is angle dependent on the lattice. The 
question of angle dependence of (J5J will be discussed else- 
where. The implied force is 



kTl 
~~2r' 



(9) 



Notice that the law itself does not depend on the stiff- 
ness, or on surface tension, but the range of validity does. 
Techniques like those in Ref. [H might be used to make 
the lower bound on r as small as possible. The reader 
might feel that the pinning mechanism is somewhat con- 
trived. If the interface were localized both at the bound- 
aries and at both the interior points, perhaps a more in- 
tuitive restriction than "floating" at the same level, then 
Z x /Z ^ a,s N 00, rendering the definition useless. 
In experiments one always has a finite system, which may 
introduce an additive term in eq:7, which is independent 
of n, but which may well diverge with N, as would be 
the case with the alternative pinning just described. But 
for any finite N, this would not come up in the force as 
in eq:7b. There are many variants here, which will be 
developed in a longer paper. 

We turn now to d = 3 and consider the van Beijeren 
BCSOS 6] model, where height changes are associated 
with vertex configurations as shown in Fig. [51 The fact 
that the allowed configurations are those of the 6— vertex 
model with two arrows in and two arrows out, guaran- 
tees height conservation around any close loop, which is 
required for the integrity of the sheet. Along the (1,0) 
direction, which is normal to the transfer direction, the 
height difference of abutting plaquettes is strictly ±1. 
Labelling the row of vertical arrows by eigenvalues of 
2Sj for spin -1/2, the following XXZ Heisenberg - Ising 
Hamiltonian 



M 



(10) 
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Sjm+i = S",a = x,y,z (cyclic boundary conditions) 
commutes with the transfer operator, provided A = 
1 - (l/2)e A '. Thus -oo < A < 1/2. The region of 
interest here is for —1 < A < 1/2 where the interface is 
rough. The ground state of (|10[) furnishes the maximum 
eigenvalue of the transfer matrix. Proceeding in the (1, 0) 
direction, the height difference at a separation r = 2nao 
for integer n is 



6h = 2j2S* 

3 = 1 



(11) 



In this letter, we are interested in the situation when 
5h = 0. The ratio of partition function with this height 
fixing to that without it is 

r*m*n) = (12) 
Z Jo ^ 

where f ' x (n) is the associated incremental free energy 
and 



f(0,2n) 



(exp 




(13) 



We are interested in the large-n asymptotics of f(8, 2n). 
We note that exp(— /3/ x (n)) is simply the probability 
that the z-component of the spin on the interval [l,2n] 
is zero. In this sense this quantity is related to Kore- 
pin's "emptiness formation probability" and its general- 
izations [2l],[22|]. To calculate exp(— /3/ x (n)) we employ 
a bosonization analysis. 
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FIG. 2: The arrow configurations of allowed configurations of 
the 6- vertex model. The rule for associating plaquette heights 
on the dual lattice is that, pointing along the arrow, the pla- 
quette on the right is higher. Vertices 1 to 4 have a diagonal 
height change in one or the other direction. Vertices 5 and 6 
have a ripple (no diagonal height change) and are deemed to 
be flat. The weight assignment, with K > 0, favors rippling 
over diagonal height change. 

It is well known that the large distance behaviour 
of correlation functions in the Heisenberg-Ising model 
HXXZ is described by a Gaussian model [23j 
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Here the scalar field $ is compact $(ac) = $(a;)+87ra, and 
the velocity v and compactification radius 4a are given in 
terms of the parameters of the lattice Hamiltonian HXXZ 



and the lattice spacing qq as a — (arccos(— A)/4ir)^' 2 ' , 
and v — ( Ja sin(47ra 2 ))/(2(l — 4a 2 )). The lattice spin 
operator Sj is expressed in terms of the bosonic field as 
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where x = jao and A is a known constant [24j . We now 
turn to the bosonization of exp^i#^™ =1 S" + ^j. We 
first note that this operator is a 2-7r-periodic function of 
which allows us to write 



Mn) = 
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Application of bosoSz then gives (the staggered part of 
does not contribute to the sum as $(x) is a slowly 
varying field) 



— / dx d x $ 
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= A(0)exp — [$(na ) - *(0)] .(17) 
\87ra / 

Here A(9) is a normalization constant. The expectation 
value of vertexop is easily calculated 



e »E?=i^ „A(9) (nc)^)", (18) 
H 



where c 1 ao is a short-distance cutoff. Using this result 
we obtain the following expression for f(6, n) 
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We can check result by comparing it to the known expres- 
sion at the free fermion point A = 0. Here f(0, n) can 
be expressed as a Toeplitz determinant, see [19l. |25| and 
references therein. The large-n asymptotics have been 
calculated in [26j with the result 
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where 



4tt 2 sinh 2 (6»V27r) 



sinh (t) 



(21) 



4 



This agrees with the general expression result. It is now 
a simple matter to calculate exp(/ x (2n)) , as for large n 
the 9- integral is dominated by its saddle point. Hence 



-0f X (n) 



A 



(22) 



where A and d are (A-dependent) constants. Thus, the 
resulting force is attractive and behaves for large r as 



kT 



(23) 



2r log(rc'/oo) 

Notice the difference between ef2 and eq:7b. 

In this Letter, we have presented exact calculations of 
the incremental free energy resulting from pinning the 
interface between two coexisting phases at two points 
separated by a distance r. The case for d = 2 is more 
detailed because it allows the interface to have a diffu- 
sive structure at the molecular level, as is known to be 
the case in real systems. The resulting force is attrac- 
tive with a remarkably long-ranged decay as 1/r. In 
d = 3, we have treated this pinning phenomenon within 
the BCSOS model of a random surface, using van Bei- 
jeren's isomorphism with the 6-vertex model. Thus, the 
interface is locally sharp. It is defined at points on a 
Z 2 square lattice. The attractive force in this case does 
not differ strongly from the d — 2 one, in that it decays 
as l/(rlogr) instead of as 1/r, a rather striking result. 
Note that these forces are of much longer range than the 
electrostatic ones induced by photon fluctuations, termed 
van der Waals 27} or the fluctuation-induced Casimir- 
type repulsion between extended objects in similar sys- 
tems. This l/(rfogr) result has also been obtained by 
Lehle, Oettel and Dietrich [lil ] in a continuum model for 
a particular, but very physically reasonable case of con- 
tact of the interface with two extended objects; these we 
may think of as spheres each with two hemispherical sur- 
faces having different wetting properties which keep the 
sphere at the interface at an orientation and immersion 
which makes the equatorial line separating the two wet- 
ting region sit exactly in the interface, as mentioned at 
the beginning of this Letter. The case we have consid- 
ered is an extreme one in that the colloidal sphere has 
been treated as though it were a point (or a unit cube on 
the dual lattice). In mitigation, our results for d=3 are 
obtained for mesoscopic modelling which retains, unlike 
the usual capillary wave scenarios, molecular scale dis- 
creteness and which also shows a roughening transition. 
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